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ABSTRACT 

Non-Maxwell i an electron velocity distribution functions, computed 
for Dupree's 1972 model of the solar transition region in a previous 
paper, are used to calculate ionization rates for ions of carbon, nitrogen, 
and oxygen. Ionization equilibrium populations for these ions are then 
computed and compared with similar calculations assuming Maxwellian 
distribution functions for the electrons. The results show that the ion 
populations change (compared to the values computed with a Maxwellian) in 
some cases by several orders of magnitude depending on the ion and its 
temperature of formation. 
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1. INTRODUCTION 

The ionization equilibrium calculations found in the literature, to 
date, make use of various theoretical and semi -empirical cross-sections 
to compute ionization and recombination rates (Allen and Dupree, 1969; 
Jordan, 1969; Summers, 1974; Jacobs et al_. , 1977a and b). All of these 
studies, however, assume that the velocity distribution function for 
ionizing electrons is a Maxwellian. 

If the ionization energy, x> is Targe compared to the local electron 
thermal energy, then it is important that we have an accurate description 

t 

of the tail of the electron distribution function in order to compute 
.the ionization rate, q^. In Table I we list the ionization energies 
for various ions found in the transition region, along with the tempera- 
ture (T ) at which the ion has its maximum population relative to n c - 
max t 

based on Jordan's (1969) results - and the ionization energy normalized 

to kT . We find that x/kT ranges from 5-15. Roussel -Dupre (1980, 
max max 

hereafter Paper I) showed that the electron distribution function is 
non-Maxwel 1 i an for energies around 6 kT, where T = 10 K, and we may 
conclude therefore that present ionization equilibrium calculations 
may need serious revision, especially for ions formed primarily in the 
region from 10 5 - 10 6 K. We also note that an increase in the popula- 
tion of the tail of the distribution function will increase the ionization 
rates, since a threshold energy is involved. On the other hand, the 
cross-sections for recombination decreases with energy and no threshold 
is required; therefore, we do not expect recombination rates to be 
affected. 

In this paper, we use the electron distribution functions, computed 
in Paper I for Dupree's (1972) model of the transition region, to compute 
ionization rates for ions of carbon, nitrogen and oxygen. We then 
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recompute the ionization equilibrium populations for these various ions 
and compare our results with calculations for a Maxwellian distribution. 


2. ' CORRECTED IONIZATION RATES 


For our purposes, we chose to work with the semi-empirical cross- 
section, given by Allen (1973), for direct ionization from the ground 
state. This cross-section is given by. 


14 .6(n,£) 
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where 6(n,£) is the number of- bound electrons with quantum numbers n and l, 
^ v (n,.£) is_ the- ionization potential .in electron volts and e is the 
electron energy. To compute the auto-ionization rates, we used the 
cross-section for excitation from state i to state j,, given by. 
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where f . . is the oscillator strength and g is a Gaunt factor 0.2). 

1 0 

These cross-sections agree with experimental measurements to within a . 
factor of 2 (See Elwert 1952). Although there are more accurate experi- 
mental results available for specific ions, we are only interested in 
estimating the effect of the electron tail on ionization equilibrium 
calculations for the ionization stages of .several elements, and approxi- 
mate' cross-sections are accurate enough for this purpose. 

In order to evaluate the effect of the tail of the electron distri- 
bution function on ionization rates, we computed the ratio of these 
rates, calculated with the non-Maxwell i an distribution function given 
by equation (22) of Paper I, to the corresponding rates calculated 
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assuming a Maxwellian distribution function. This ratio, for direct 
ionization from the ground state is given as a function of I(= x/kT) and 
u , the critical speed above which deviations from Maxwellian are large 

V * 

by. 
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For auto-ionization we find. 
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where e^c(x)is the complementary error function; E^ Q {x) and E^(x) are 

p 

the tenth and first exponential integrals respectively and G(I,u ) is 
defined by. 


6(1, u^) - 


-J 


u 4 


exp(u^x) • 


dx . 


By simple substitution, we can. show that equations (3) - (6) satisfy the 
conditions. 

Rf (I» u c ) 1 as 


R^(I.I) = R- (IJ) 

R a ( 1 * u c ) "*■ 1 as u . c ■*" “ 


R a (I = R a (IJ) * 

As a typical example, we list, in Table II, values of and R Q 
for the carbon IV ion, -as a function of ' temperature and the critical 
speeds, u c , computed from Dupree's model in Paper I. The non-Maxwell i an 
tail leads to ionization rates which can be as much as a factor of ten 
larger than those computed assuming a Maxwellian. In addition, since 
radio observations predict a lower pressure (or density) by a factor of 
three or more for the transition region than 'EUV observations (cf., Dulk 
et aK, 1977), we have also listed these ratios for Dupree's model with 
the pressure decreased by a factor of three. The latter results yield 
larger ratios simply because the critical speeds are decreased, extending 
the tail over a larger range in velocity space and further increasing 
the population of high energy electrons relative to a Maxwellian. 
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3. IONIZATION EQUILIBRIUM CALCULATIONS 
We are interested in computing the effect of these higher ionization 

rates on the relative populations of the various ionization stages of a 

* 

given element in ionization equilibrium in the transition region. The 
relative population of two neighboring ionization stages for a given 
element is computed .directly from 


n(x +l,tf1 ) _ j < ’i (x ' fr ") 
n(x +m ) E r c. r (x +m ) 


(7) 


where q^(x ) represents the ionization rate for an ion in ionization 
stage m by a process i, a^(x +m ) represents the recombination rate for an 
ion x +m by a process r and the sums i and r are over all ionization-and 
recombination processes, respectively. 

In our analysis, we considered direct ionization from the ground 
state, auto-ionization, radiative recombination and di-electronic 
recombination to be the most important ionization and recombination 
processes for ions in the transition region. The di-electronic recombination 
rates used in our computations were used in the calculations of Allen and Dupree 
(1969) and provided by Joselyn. The states which auto-ionize for a given 
ion were obtained by means' of a scheme described by Jordan (1969). The 
ionization rate from the ground state was computed in the manner described 
in the previous section and the radiative recombination rates were 
computed' directly from equation (2-13) of Allen and Dupree (1969), assuming 
a Maxwellian distribution function; however, this is a good approximation 
since the tail only affects recombination rates by a few percent at 
most. 

Finally, we computed the .fractional ion populations (N^/N g ) for the 
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various ionization stages of carbon , nitrogen and oxygen in ionization 
equilibrium. The results, using Dupree's model to compute u c , are 
presented in Figures 1-3; while the results, using Dupree's model with 
the pressure reduced by a factor of.thr.ee to compute u c , are shown in 
Figures 4-6. In each of these- figures, the solid curves are based on 
the assumption that the distribution functions are Maxwellian. The 
dashed curves take into account the effect of the electron tail. In the 
region where the temperature gradient is steepest, we find that the 
lower ionization stages of all these elements are. depleted by the 
increased ionization caused by the tail. The higher ionization stages 
have increased populations throughout the upper half of the transition 
region suggesting that ions, thought to form only in the corona, can 
have substantial populations in the transition region. The changes in 
population can be many orders of magnitude for Dupree's model and are 
even higher when the pressure is lowered by a factor of three. 

As pointed out in Paper I, it is impossible to pinpoint the exact 
value of the critical energy. As a result, we decided to present 
calculations for a minimum and maximum critical energy. The results 
presented above correspond to the minimum value. In Figure 7, we present 
ionization -equilibrium populations for carbon in the same format described 
above, computed with the maximum critical energy. We find that these 
results show no apparent difference from those obtained with a Maxwellian 
for the majority of ionization stages of carbon except C VI which has a 
large ionization potential. 

We can conclude that the effect of the high energy tails on ionization 
equilibrium is sensitive to the location of the tail. This means that 
the exact magnitude of this effect cannot be pinpointed by our method of 
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approximation and that this problem will remain unresolved until detailed 
calculations, which would involve a solution of the Boltzmann-Fokker- 
Planck equation, are presented. However, we have argued in Paper I that 
the minimum critical energies used' in our calculations represent a good 
approximation for the location of the tail. - 

• , 4. SUMMARY 

The effect of the high energy tails on ionization equilibrium 
calculations is summarized in Figures 1 - 7. We found that the ion 
populations change (compared to the values computed with a Maxwellian) 
by as much as several . orders of magnitude depending on the ion and its 
temperature of formation. We also found, however, that these results, 
are sensitive to the location of the high energy tail in velocity space 
and that detailed solutions of the Boltzmann-Fokker-Planck equation are 
necessary to determine the exact magnitude of these effects. However, 
arguments presented in Paper I suggest that the minimum critical energies 
represent a good approximation for the locations of the tails so that 
the results obtained with these critical energies are reasonable approximations 
for the changes -in ion populations. 

The effect of the high energy tails on ionization equilibrium 
calculations is significant. This result has serious implications for 
the interpretation of EUV line profiles to yield models for the transition 
region. The high energy tails affect the collisional excitation rates 
(not computed here), ion populations and the temperature range over 
which ions form. Both excitation rates and ion populations' enter into 
line intensity calculations. We can conclude that our results may 
significantly change the models which we compute from EUV line observations. 
However, since detailed calculations which include the effects of high 
energy tails on computed line profiles have not been presented, it is 
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difficult to conclude, at this point, exactly what changes will result 
in present models. Calculation of a self-consistent model, for the 
transition region, which reproduces EUV line observations is a topic for 
future research. 

Finally, we note that the effect of magnetic fields at an angle to 
the temperature gradient was not included in our calculations. Consider- 
ing the multitude of magnetic field configurations observed on the solar 
surface, a study of transport processes in the presence- of magnetic 
fields and large temperature gradients is essential and remains a 
subject for future research. In addition, we mention that the analyses 
presented here and in Paper I are not limited to the transition region 
but are applicable to any plasma where steep gradients are encountered. 
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TABLE I 

RATIOS OF IONIZATION ENERGIES' TO LOCAL TEMPERATURES OF 
FORMATION FOR SEVERAL IONS 


Ion 

Ionization Energy 

Log 

x/kT 



X (ev) 

3 max 

max 

C 

ii 

24.4 

4.30 

1.42(1} 

C 

in 

47.9 

4.85 

7.84 

n 

ri 

29.6 

4.40 

■ 1.37(1) 

n 

hi 

47.5 

4.95 

6.17 

n 

V 

98.9 

5.30 

5.75 

0 

in 

54.9 

4.90 

8.01 

0 

IV 

77.4 

5.25 

.5.05 

0 

V 

114.0 

5.40 

5.26 

0 

VI 

138.0 

5.50 

5.06 
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TABLE II 

RATIOS OF CIV IONIZATION RATES COMPUTED WITH A TAIL 
' TO THOSE COMPUTED WITH A MAXWELLIAN 


Log T 

R i/R i * 
P = P o 

R a/R * 

a 

p .p o 

R i/R.j* 

P = P o /3 

V.* 

P * P 0 /3 

4.7 

1.01 

1.00 

3.32 

1.58(1) 

4.8 

3.13 

1.26(1) 

3.03(1) 

3.08(T) 

5.0 

3.64 

8.87 

1.09(1) 

1.39(1} 

5.2 

1.80 

5.42 

2.99 

6.02 

5.4 

1.06 

1.01 

1.31 

1.05 

5.6 

1.01 

1.00 

1.04 

1 .01 

5.8 

1.00 

1.00 

1.01 

1.00 

6.0 

1.00 

1.00 

1.00 

1.00 

6.2 

1.00 

1.00 

1.00 

1.00 
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Figures 1-3 


Figures 4-6 


Figure 7 


FIGURE CAPTIONS 

Fractional ion populations for the ionization 
stages of carbon, nitrogen and oxygen (Dupree's 
model). 

Fractional ion populations for the ionization 
stages of carbon, nitrogen and oxygen (Dupree's 
model vn'th pressure reduced by a factor of 
three). 

Fractional ion populations for the ionization 
stages of carbon (U^ - (0.4)^(U^) D ). 
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